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Abstract 

Classical Stokes' drift is the small time-averaged drift velocity of suspended non- 
diffusing particles in a fluid due to the presence of a wave. We consider the effect 
of adding diffusion to the motion of the particles, and show in particular that a non- 
zero time-averaged drift velocity exists in general even when the classical Stokes' drift 
is zero. Our results are obtained from a general procedure for calculating ensemble- 
averaged Lagrangian mean velocities for motion that is close to Brownian, and are 
verified by numerical simulations in the case of sinusoidal forcing. 

PACS numbers: 02.50-r, 05.40+j, 05.60+w 

A travelling wave in a fluid gives suspended particles a small drift velocity known as 
Stokes' drift |1|, [2|, |3[]. When there is more than one wave, the drift velocity is calculated by 
summing the contributions from each wave [|], ^ . In this letter we consider the influence of 
diffusion on the magnitude and direction of the drift velocity. As in the classical (diffusion- 
less) case, the amplitude of the travelling wave is assumed small compared to its wavelength; 



a non-zero drift velocity appears at second order in the amplitude. In the presence of more 
than one wave, the classical Stokes' drift can sum to zero. Diffusion then produces a non- 
zero drift velocity whose magnitude and direction depends on the diffusivity of the suspended 
particles. 

Several mechanisms for the directed motion of small particles without a net macroscopic 
force have been proposed in the last 10 years |6|, |7|, §, || [L(J. Interest in such 'ratchet' 
effects has been motivated by the search for the mechanisms of biological motors, such 
as the conversion of chemical energy into directed motion by protein molecules, and by 
possible applications, such as the separation of particles in solution based on their diffusion 
coefficients. In both these cases small particles are believed to follow dynamics that are 
overdamped (first derivative in time) and noise-dominated. A drift velocity dependent on the 
size of suspended particles in solution has been produced experimentally using an asymmetric 
periodic potential turned on and off periodically ||11|| . Published theoretical models [12, [13 
combine a periodic asymmetric potential in one dimension with non-white fluctuations. 

In this letter we consider motion in arbitrary dimensions that is diffusion-dominated. 
There is also a small deterministic forcing whose amplitude will be used as an expansion 
parameter; a drift velocity appears at second order and depends on the diffusivity. Thus 
diffusion due to microscopic motions, for example diffusion of particles in solution, can be 
exploited using a carefully-chosen combination of forcings to produce a net motion that 
depends on the diffusivity. We illustrate the effect with sinusoidal forcing and compare our 
calculations with numerical results in one and two space dimensions. It is possible to arrange 
the wave motions so that particles of different diffusivities have a time-averaged drift velocity 
in different directions, resulting in what we call 'fan-out'. This may have applications for 
sorting particles according to their diffusivities. We show numerically that the fan-out can 
have an angular range of more than 180 degrees. 

We first develop an expansion scheme for motion that is overdamped and diffusion- 
dominated. Consider a stochastic process X = (Xf) t > taking values in IR m and satisfying 
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the following stochastic differential equation [TJ], KB : 



dX t = efC& t ,t)dt + dW t , 0<e<l. (1) 

The vector X t is the particle position at time t. Its ensemble average, to be denoted below 
by angled brackets, is the Lagrangian mean position at time t. W is an m-dimensional 
Brownian motion, with W = and (W 4 ■ W 4 ) = ma 2 t, i.e. W represents a purely diffusive 
motion, with diffusivity 

D = \a 2 . (2) 

The remaining term in (|I|) is the deterministic forcing, a function of Eulerian position x and 
time t: 

/:l m xR + ^ R m . (3) 

The real constant e satisfies < e < 1. 
We now expand in powers of e. Let 

X t = X| 0) + eXf } + e 2 x| 2) + . . . , with initial condition X = 0. (4) 

The leading terms of the stochastic equation of motion, equation (Eh, are as follows. 



e°: 



e 1 : 



dxj 0) = dW t> (5) 

giving 

Xi 0) = W t . (6) 
dX| 1) = /(x| 0) ,t)dt, (7) 

giving 



X| X) = [ f(W s ,s)ds. 
Jo 
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dXf ) = (X t (1) .V)/(x| 0) ,t)dt, (9) 
giving the second-order drift velocity as 

A x < 2) = jT(/(W., a) ■ V)/(W t , f) da, (10) 

where V/ is the spatial gradient of /. 

In the classical derivation of Stokes' drift there is no motion at zeroth order [jlj [|. Here 
the motion at zeroth order is purely diffusive, with (X.^) = f° r an y positive time. When 
f(x,t) is a sum of functions that are periodic in t at any fixed x, the drift velocity also 
vanishes at first order in e. That is 

lim i(xl 1} ) = 0. (11) 

At second-order the drift velocity is in general non-zero. It is given by the following ensemble 
average: 

U = e 2 lim i<XtW) (12) 

t^co t 

In one space dimension the expression ( |i~3"|) reduces to 

(14) 
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where f'{x,t) = £f(x,t). 

We now consider the case where the deterministic forcing is a sum of sinusoids: 

n 

f(x,t) = A e k e cos(k e ■ x - u e t), (15) 

where Ag and ui are constants. The vector ki defines the direction of propagation of wave I. 
The drift velocity flT3| ) for this case is 

n „t 

U = e 2 Y,\^A]\h\ 2 h (sin(-^-(W t -W s )+o;^-s))))d S 



n rt 

^e 2 Y^ lim Aj\ki\ 2 ki / sm(u)i(t — s))e~ 



\ki\ 2 D(t-s) di 



(16) 
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Each wave makes a contribution to the drift velocity in its direction of propagation. For 
D — 0, the weighting factor is proportional to the square of the amplitude. This is the 
classical result obtained by a transformation from Eulerian to Lagrangian coordinates [0], . 
In the case of a surface wave over deep water, the first order motion of a suspended particle is 
a circle with radius A(k^/ui^ the quantity A\ \k^\ 2 k^/ui is proportional to the time-averaged 
momentum per unit area || In the presence of diffusion, the contribution from wave t 
is reduced by the dimensionless factor (1 + af) , where ote = D\k(\ 2 uJ l . Diffusion reduces 
the Stokes drift due to any one wave by smearing out the distribution of particles, working 
against the tendency of particles to spend longer in regions where the force acts in the 
direction of propagation than in those where the force acts in the opposite direction. The 
attenuation is strongest for waves with large wavenumbers or small velocities. 

Dependence of drift velocity on diffusion can be exploited as follows: there is in general 
a non-zero drift velocity due to diffusion even when the classical Stokes' drift is zero. We 
write the drift velocity (|1(J) as a sum of the classical Stokes' drift and a diffusion-dependent 
contribution: 

U = U + U s , (17) 

where 

U = U\ D=0 = ^±A 2 \k £ \ 2 ^ (18) 
i=i u t 

and 



A]\h\ 



2 k e aj 



(19) 



uji 1 + aj 

The classical Stokes drift Uq can be made to vanish by choosing a forcing f(x, t) consisting 
of two wave trains propagating in opposite directions. For the latter example, we can work 
in one space dimension, defined as the direction of propagation of wave I = 1: 

f(x, t) = A x ki cos(fciX — u>it + 0i) + A 2 k 2 cos(k 2 x — u 2 t + 2 )> (20) 



where A iy ki, u>i and (pi (i = 1, 2) are constants and kik% < 0. For simplicity, we suppose 
that k± 7^ ±^2 and u)\ ^ ±u>2~, this avoids cross-terms in the classical Stokes' drift. The drift 
velocity including diffusion is then given by 

u= ^(Afi( 1+al y +Al k i (1+a f } -y (21) 

To set U — requires Afk^/ui = —A\k\j 'u; 2 . Then U = U s where 

C/s^le^ffl + a?)" 1 - (l + ai) _1 V (22) 



For large diffusivity the drift velocity tends to zero because the contribution of each wave 
tends to zero. Thus there is an intermediate value of diffusivity that maximises U s . If the 
forcing frequencies and wavenumbers are fixed and Uq = 0, the drift attains its maximum at 
the value of D satisfying a±a2 = 1. 

Figure [j] shows the drift velocity as a function of diffusivity with the forcing a sum of two 
sinusoids for a choice of parameters that gives Uq = 0. In Figure ^] the calculated drift is 
compared with numerical results, with the same choice of parameters and D = 0.125. The 
solid line in Figure |](a) is the mean value of X t as a function of time, averaged over 10000 
numerical realizations of the stochastic differential equation ([I]), and the dotted line is Ut, 
with U given by (E^j). In Figure ^(b) we show, as a function of time, the difference between 
the numerically-calculated mean displacement and Ut. Figure |2|(c) demonstrates that the 
motion is close to Brownian; a histogram of values of X 4 at t — 1000, R(y), is compared 
with the Gaussian probability density function with mean Ut and variance a 2 t (solid line). 

In general the expressions ( |16D and (|19D are vector relations. Thus, in more than one 



space dimension, the direction as well as the magnitude of the drift velocity depends on the 
diffusivity, producing fan-out of the drift velocity vectors. We illustrate this effect in Figure 
H|, constructed with the forcing being a sum of four sinusoids in two dimensions. In (a), the 
vector Aikt is shown for each of the four waves. The parameters are A\ = 1.0, A2 = 0.8, 
A 3 = 0.7, A A = 0.7; k x = (1.0,0.0), k 2 = (2.0,-4.0), k 3 = (-3.0,0.7), k 4 = (-0.96,4.56). 
We take uo = vk with v = 1. Figure |^(b) depicts the fan-out in the directions and magnitudes 
of the drift velocities for nine different values of diffusivity. Each arrow is U s for one of the 
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following values of D: D = 0.1 (leftmost arrow), 0.2... 0.9 (rightmost arrow). For larger 
values of D, the direction of U approaches more closely that of k\. 

The fan-out effect shown in Figure |3] is due to the different rates at which the contribution 
from waves decreases as the diffusivity is increased, destroying the exact cancellation imposed 
at D = 0. More light is shed by considering the small-diffusivity and large- diffusivity limits 
of ©. 



1. If D\k e \ 2 /uj e «1 VI then 



\h I 6 

A 2 t ^-k E + 



(23) 



2. If D\k e \ 2 /uj e > 1 V/ then 
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(24) 



In the limit of small diffusivity (|l|) the drift velocity is proportional to D 2 and the direction 
is approximately opposite to that of the wave with the largest value of y4 2 |/c| 6 c<j _3 . In the 
opposite limit (0) the drift velocity is proportional to D~ 2 and the direction is approximately 
parallel to that of the wave with the largest value of y4 2 tu|A;| _2 . 

In summary, we derive a general expression for the drift velocity of diffusing particles from 
a stochastic asymptotic expansion scheme for motion that is Brownian plus a small deter- 
ministic forcing. The drift velocity is in general non-zero even when the classical Stokes' drift 
vanishes. For example, several counterpropagating sinusoidal forcings produce a drift veloc- 
ity that depends on the diffusion coefficient and the intensities, frequencies and wavenumbers 
of the forcings. Thus, given a collection of particles with different diffusivities, the deter- 
ministic forcings can be tuned to separate particles of a particular type by optimizing their 
stochastic Stokes' drift. 
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Figure 1: Stochastic Stokes' drift in one dimension. There is a non-zero drift velocity due 
to diffusivity even though the classical Stokes' drift vanishes, due to the different rates at 
which the contributions from each wave decrease as the diffusivity is increased. The drift 
velocity, fl22[) , is given as a function of diffusivity for e = 0.1, Ai = k\ = uj\ — 1, k 2 = —2.42, 
u 2 = 0.47. 
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Figure 2: Stochastic Stokes' drift: comparison with numerical solution for sinusoidal forcing 
in one dimension, (a) Mean value of X t . (b) Difference between the numerical mean value 



and the second-order result (|22|) . (c) Distribution of X t at t = 1000, a = 0.5. 
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Figure 3: Stochastic Stokes' drift in two dimensions, (a) Directions and magnitudes of 
the four sinusoidal forcings. The classical Stokes' drift (p~8|) is zero, (b) Resulting stochastic 
Stokes' drift as a function of diffusivity. Here D = 0.1, . . . , 0.9 (largest D on the right). In 
the axis labels, the subscripts indicate vector components. Note the fan-out of more than 
180 degrees. 
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